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I. INTRODUCTION 

Finding exact solutions of the metric field equations 
of General Relativity (GR) is, in general, a non-trivial 
task, though over the years a wide collection of such so¬ 
lutions has been obtained [l[ . Among them, we find those 
representing the collapse of spherically symmetric bodies 
(Schwarzschild) with charge (Reissner-Nordstrom) and 
rotation (Kerr-Newman). These solutions put forward 
the dramatic effects that massive bodies and electric 
charges might have on the causal structure of spacetime. 
The Schwarzschild black hole possesses a null hypersur¬ 
face, called event horizon, that prevents any particle that 
goes through it from coming out again. The effect of 
electric charge is also remarkable, since a second event 
horizon may appear inside the black hole giving rise to a 
much more complex causal structure, as is manifest from 
the corresponding Penrose diagram. 

Another relevant aspect of the internal structure of 
black holes resulting from gravitational collapse is the 
existence of a singularity at their center. Within GR 
this is an unavoidable consequence provided that i) there 
exists a (future) trapped surface, ii) the matter energy- 
momentum tensor satisfies reasonable energy conditions, 
namely, the null energy conditioifl, and iii) global hy- 
perbolicity holds 0 (see also i). Though the very 
definition of the concept of spacetime singularities re¬ 
mains elusive, it is traditionally seen as related with the 
existence of a troublesome region of spacetime marked 
by the divergence of some geometric magnitudes. A 
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^ Note that in the original formulation of the theorems this is 
rather a geometric statement on the positivity of the Ricci tensor 
for null vectors N°‘, namely, RabN°‘N^ > 0, which, in the case 
of GR via the Einstein equations, becomes a statement on the 
energy conditions. 


standard way to characterize them is to consider in¬ 
variant polynomials constructed from the Riemann ten¬ 
sor R^vap and see if they blow up somewhere. In 
this sense, in the Schwarzschild case, the Kretschmann 
scalar K = becomes Ks = whereas 

in the Reissner-Nordstrom (RN) case we have AT/jat = 

, where rg = 2M is the Schwarzschild 
radius and = 2Gq^ is a length scale associated to the 
charge. The higher degree of divergence as r —^ 0 in 
the charged case suggests that the energy associated to 
the electric field contributes to worsen the pathological 
behavior of the geometry as the sources are approached. 

A more powerful characterization of spacetimes con¬ 
taining singularities, however, is provided by the notion 
of geodesic completeness, namely, whether an affine pa¬ 
rameter on every geodesic curve extends to arbitrarily 
large values or not. In a singular spacetime, there ex¬ 
ist geodesic curves for which the affine parameter cannot 
be extended to arbitrarily large values, i.e., they start or 
terminate at a finite value of the affine parameter. This 
captures the intuitive idea of a spacetime singularity as 
the phenomenon by which “an observer’s future comes 
suddenly to an end”. If the affine parameter can take 
arbitrarily large values on the real line, then we say that 
the spacetime is geodesically complete and nonsingular 
regardless of the behavior of the curvature invariants [J . 
This is generally accepted as the most reliable criterion 
to determine if a spacetime has a singularity d, d , and 
will be the one accepted in this work. 

To achieve singularity avoidance, one is thus forced 
to remove at least one of the assumptions on which the 
singularity theorems are based. For example, hypotheti¬ 
cal matter-energy sources violating the null energy con¬ 
ditions have been introduced in the literature. Among 
these we find phantom black holes 0 or solutions sup¬ 
ported by non-linear electrodynamics Q, which are able 
in some cases to obtain regular black hole solutions. In 
this work, however, we shall take another way round and 
accept the viewpoint that the divergent behavior of cur- 
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vature scalars as r —>■ 0 is an indication that the classical 
description provided by GR breaks down in that region, 
and that some improved theory of the gravitational field 
should be used to correctly describe the physics in the 
innermost regions of black holes. In fact, it is widely as¬ 
sumed that the quantum effects of the gravitational 
field should manifest themselves at curvature scales of 
order K ~ 1/^p, where /p = hG/c^ is the Planck length 
squared, and thus modify the classical description pro¬ 
vided by GR. Given our current limited understanding 
on quantum gravity, its impact on the geometry around 
black hole singularities is difficult to foresee, though some 
relevant results have been obtained recently using pow¬ 
erful quantization techniques in simplified scenarios [HI- 
A different approach to this fundamental question con¬ 
sists on assuming that the bulk of the quantum effects 
of gravity can be captured by some effective theory in 
which the undesirable features of spacetime singularities 
can be avoided [l^. In fact, since a basic identifying fea¬ 
ture of quantum gravity is to cure spacetime singularities, 
it seems reasonable to expect that some effective classi¬ 
cal geometry without the shortcomings of GR should be 
recoverable in some low energy regime. It is in this phe¬ 
nomenological context that this paper is framed. We note 
that other approaches to this problem exist in the litera¬ 
ture such as those inspired on non-commutative quantum 
gravity [l^, variations of Newton’s constant such as in 
RG gravity [l^, modifications of the dispersion relation 
+ rn? within the so-called Gravity’s rainbow 
(see [13 for the original proposal and [l^ for black hole 
solutions), the non-Lorentz invariant (at high-energies) 
Horava-Lifshitz gravity , etc. 

In a number of previous works [isjl , it has been shown 
by some of us that the innermost structure of spher¬ 
ically symmetric, electrically charged systems coupled 
to certain metric-affine extensions of Einstein’s theory 
might be completely free of curvature divergences. This 
happens, in particular, when quadratic curvature terms 
are added to the Lagrangian, and also for Born-Infeld 
type modifications of the gravitational theory [I^. The 
metric-affine (or Palatini) version of those theories is gov¬ 
erned by second-order equations and is ghost-free. This 
allows to obtain exact analytical solutions and prevents 
severe shortcomings that arise in the standard metric 
(or Riemannian) formulation of those theories. The Rie- 
mannian approach assumes that the connection should 
be metric-compatible (and thus given by the Ghristoffel 
symbols of the metric), which generically leads to higher- 
order equations for the metric. In the Palatini approach 
this a priori constraint is relaxed and the connection 
is determined through the field equations (see e.g. [ 2 ^ 
for a review on Palatini gravity), which keeps the equa¬ 
tions second ordei0. The resulting equations are more 
tractable and allow to find exact solutions in some cases 


of interest and explore their physical properties without 
resorting to perturbative treatments. One then finds 
that, for configurations extending the RN solution to 
this framework, the central curvature divergence can be 
avoided when a certain charge-to-mass ratio is satisfied. 
For arbitrary values of the charge and mass parameters, 
however, curvature divergences arise on a sphere of area 
A = 47rr^, with = l^r^, being the length scale char¬ 
acterizing the high-curvature corrections in the gravity 
Lagrangian [see (1101) below]. These divergences are much 
milder than in the GR case, dropping from K r^ /r® to 
K ^ l/(r — Tc)^ as r —)> rc is approached. This change 
occurs in a smooth but non-perturbative way. It turns 
out that the surface r = rc represents the throat of a 
wormhole, a tunnel to another region of spacetime, gen¬ 
erated by the interplay between the electric field and the 
Palatini gravity. As we will see, the existence of this 
wormhole, a topologically non-trivial structure, has deep 
implications for the understanding of the properties of 
curvature divergences and their relation with spacetime 
singularities. 

The aim of this work is to progress in the understand¬ 
ing of the geometric properties of these solutions and 
the physical meaning/implications of the existence or ab¬ 
sence of curvature divergences. Motivated by the fact 
that smooth solutions, which are geodesically complete, 
can exist arbitrarily close to solutions with curvature di¬ 
vergences, we explore in detail the geodesic structure of 
these spacetimes for the whole space of configurations of 
mass and charge. One of our goals is to determine if 
the existence of curvature divergences at r = rc implies 
that geodesic curves terminate there in a finite affine pa¬ 
rameter (geodesic incompleteness). We hnd the answer 
to this question to be negative, i.e., that despite hav¬ 
ing curvature divergences at r = rc geodesics can be 
smoothly extended through that region. This fact puts 
forward that a spacetime can be non-singular (geodesi¬ 
cally complete) despite having localized curvature diver¬ 
gences, which calls for a reconsideration of the role typ¬ 
ically attributed in the literature to curvature invariants 
for the characterization of spacetime singularities. This 
is the main result of this work. 

Our approach, therefore, consists on studying whether 
a given spacetime, specified below in Eqs.([T])-(l8]), is sin¬ 
gular or not using standard methods developed in well- 
established classical literature on the theory of General 
Relativity. In particular, we are following Geroch’s anal¬ 
ysis 0 , from which he concluded that geodesic complete¬ 
ness is the fundamental criterion above other intuitive 
criteria such as that of considering divergences of curva¬ 
ture scalars. The crucial point of that approach is that 
geodesics exist and be complete. The underlying reason 


^ Note that metric-afBne geometries seem to be of relevance for 
the proper description of solid state physics with defects on their 


microstructure such as Bravais crystals or graphene |21| . This 
yields promising new avenues for our understanding on the mi¬ 
croscopic structure of spacetime and gravitational phenomena 
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is that the existence of geodesics can be interpreted as 
the existence of physical observers. The fact that those 
observers may experience intense tidal forces or defor¬ 
mations in some regions due to curvature divergences is 
irrelevant as long as they exist. In a consistent theory, 
physical observers should always be well defined (com¬ 
plete geodesics), not disappear at some future time or 
come into existence at a given instant, which are exam¬ 
ples of incomplete geodesics. 

We note that the criterion about geodesic completeness 
is theory-independent. Whether a given spacetime is sin¬ 
gular or not is independent of where that given geometry 
originates from, i.e., it is irrelevant if it is a solution of 
Einstein’s equations or of any other physical theory. For 
this reason, we neither discuss in detail the derivation 
of the spacetime metric whose geodesics are analyzed in 
this work nor the different models that give rise to that 
solution, which have been the subject of previous works 
mentioned above. 

The paper is organized as follows: in Sec|nl we intro¬ 
duce the background geometry we are interested in, and 
describe its wormhole and horizon properties. The Eu¬ 
clidean embeddings of this geometry are constructed in 
Sec lIIIl and the conformal diagrams in Sec lIVI In SecIVl 
we provide a detailed description of the geodesic struc¬ 
ture for null and time-like curves and different configu¬ 
rations of the mass and charge parameters. We conclude 
in Sec I VII with a summary and some future perspectives. 


II. BACKGROUND GEOMETRY 


The geometry we are interested in has been derived in 
detail in a number of previous works [l^ and takes a par¬ 
ticularly simple form in ingoing Eddington-Finkelstein 
coordinates 


= —A{x)dv^ -\ - dvdx + r^{x)d^^ , 


where 


A{x) = 




1 - 


rs (I -b SiG{r)) 


1/2 



( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


being Xc a constant defined as Xc = where = 

2GNq^ is a length scale associated to the electric charge 
that, together with the Schwarzschild mass, Mq = rs/2, 
characterizes the solution. The scale 1^ characterizes 
the high-curvature corrections in the gravity Lagrangian. 


The function G{z), with z = r/cc, can be written as an 
inhnite power series expansion of the form 

G[z) = - -b [f3/4iz) + fr/iiz)] , (6) 

where fxiz) = — z^^] is a hypergeometric 

function, and Sc ~ 0.572069 is a constant. For z ^ 1, 
G(z) ~ — I/-Z yields the expected RN solution of GR, 
with (Tj- Ri I, x^(x) ~ x^, and 

rc Xg 

+ ( 7 ) 

The location of the horizons in this spacetime is almost 

coincident with the predictions of GR except for config¬ 
urations with small values of the parameters xs and Xg 
(microscopic black holes) (l^ . 

With a local redefinition of the time coordinate, dv = 
dt -\- dx/ (Acr+), the line element ([IJ can be written as 

ds^ = —A(x)dt^ + . dx^ + x^(x)dfl^ , (8) 

B{x) 


with B{x) = A(x)a‘^. Note that one could absorb the 
factor (T-|_ into a redefinition of the coordinate x to turn 
the line element into a more standard Schwarzschild-like 
form with B[x)~^dx^ = A(x)~^dx^. Such a replace¬ 
ment, though totally valid, would spoil the simple rep¬ 
resentation of x'^{x) introduced in ([5]). It must be noted 
that the coordinate x is defined on the whole real axis, 
x s] — oo, -boo[ (see Fig. [T]). As a result, one can readily 
see that the area function S = 47rr^(x) has a minimum of 
size Smin = dTrr^ at x = 0. The existence of this minimal 
sphere signals the presence of a wormhole (see e.g. 
and [ 2 ^ for references on the topic). This wormhole is 
a non-trivial topological structure supported by a sp her- 
ically symmetric electric field without sources 1^ . 
which can be interpreted as a geon in Wheeler’s sense 


On the other hand, the possibility of using the func¬ 
tion r(x) as a coordinate is subject to an important re¬ 
striction. In fact, since dx^ = a^dx^ ja-, the change of 
coordinates is ill defined at x = 0, where r = rc, be¬ 
cause dx/dx = 0 at that point. Therefore, the use of r 
as a coordinate is only valid in those intervals in which 
r(x) is a monotonic function [J. According to this, one 
would need two copies of the coordinate r(x) to cover the 
whole range of x, one for the interval in which r grows 
with growing x and another for the interval in which r 
decreases with growing x (where dx = —a+dxj^JWZ). 

Insisting on the use of r as a coordinate leads to an 
interesting effect related with the advanced Eddington- 
Finkelstein coordinate v used in (HD. For null and time¬ 
like radial geodesics, we have ds^ < 0, which implies 


„ 2 

— Adv^ -\ - dvdx < 0 . (9) 

Inside the event horizon, A < 0 implies that dx < 0, i.e., 
all null and time-like geodesics move in the decreasing di¬ 
rection of X as n moves forward in time. Now, given that 
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r{x) 



consequence of Palatini theories with higher-curvature 
corrections. 

In these theories, metric and connection must be varied 
independently in the action to obtain the field equations. 
One then finds that the connection can be solved as the 
Levi-Civita connection of an auxiliary metric related 
with the spacetime metric via the expressions 


^ ^ Vdet S) , (11) 

Vdet S ^ ^ 


where 


S = 


0'-l2x2 0 

0 ^+12x2 


( 12 ) 


Figure 1. Representation of r{x) (solid curve), defined in 
(0, as a function of the radial coordinate x in units of the 
scale Vc- The dotted lines represent the function | 2 ;|. 


r^(x) has a minimum at a: = 0, the relation between dx 

1 /2 

and dr in the cc > 0 sector is dx = dra+la_ , whereas 
in a: < 0 it is dx = —drcr+/(T_ . Therefore, ingoing 
geodesics inside the horizon, for which the evolution al¬ 
ways goes in the decreasing dx direction, propagate in the 
decreasing direction of the area function r^(x) if x > 0 
but in the growing direction if x < 0. This means that 
for an observer using r as a coordinate, ingoing geodesics 
move towards the wormhole if x > 0 but away from it 
if X < 0. For outgoing geodesics we find an analogous 
effect. 

The line element (jH]) was first obtained in as a 
solution of the combined system of Maxwell’s equations 
for a spherically symmetric, sourceless electric field plus 
a quadratic extension of GR with action 

Squad = ^ / [i? + |2(aE2 ^ 

" i (10) 

where = SttG, a is a dimensionless constant, g 
is the determinant of the spacetime metric R = 
g^'" and Rfi^iT) is the Ricci tensor of the con¬ 
nection r = which is a priori independent of the 
metric gap (metric-affine or Palatini formalism), and 
= dfj.A„ — dvA^ is the field strength tensor of the 
vector potential A^. Torsion, = (P^j, — is 

set to zero for simplicity [s^. Quadratic actions of the 
form (uni) are motivated by well established results from 
the theory of quantized fields in curved spacetimes [3l|. 
Remarkably, the line element ([8|) is also an exact solution 
of the Born-Infeld g:rav ity theory originally proposed by 
Deser and Gibbons and investigated in further detail 
in M- Wormholes with similar properties can also be 
found in the simpler scenarios of f{R) theories [s^ when 
nonlinear couplings in the electromagnetic field are con¬ 
sidered. This fact suggests that wormholes are a generic 


is a deformation induced by the energy density of the 
electric field. Here a± = 1 ± and I 2 x 2 represents 

a 2 X 2 identity matrix. From the metric field equations, 
one finds that the line element defined by h^i, takes the 
form 




—C{x)dt^ 



dx^ + x^dfl'^ . 


(13) 


Using the relations (ED, one readily verifies that G(x) = 
Aa+, with A defined in ©, and that r^ = x^ja which 
explains the dependence on x of in Eq.([5D. 

As shown above, the line element ([5]) [and also (fT^ ] re¬ 
covers the general relativistic Reissner-Nordstrom (RN) 
solution for r ^ Xc. However, as r —>■ Xc (equivalently 
X —^ 0) one finds important departures from the RN so¬ 
lution. A series expansion of the function A(x) indicates 
that depending on the value of the charge-to-mass ra¬ 
tio, di, the behavior of the solutions might differ sub¬ 
stantially. In fact, defining the number of charges as 
Nq = qje, where e is the electron charge, we have 


lim A(x) 

r-^-Vc 


Nq (di - Sc) I rc Nc - Nq 

4:Nc Side \J r - rc 2Nc 


+ O {y/r - rc) , 


(14) 


which shows that the metric is finite at x = Xc for = 5c 
but diverges otherwise. For convenience we have intro¬ 
duced the constant Nc = yj^jacm ~ 16.55, where acm 
is the fine structure constant. The smoothness of the 
geometry in those configurations with di = 5c and the 
absence of sources that generate the electric field are cru¬ 
cial elements to confirm that the coordinate x is defined 
over the whole real axis. Once this is accepted, a worm- 
hole structure arises which naturally explains the electric 
charge of the solutions as a topological effect. The re¬ 
sulting object can thus be interpreted as a geon [2^, a 
self-consistent gravitational-electromagnetic entity, with 
a non-trivial topological structure 

A careful analysis of the horizons in these geometries 
reveals the existence of the following different classes of 
solutions, according to the charge-to-mass ratio, (5i, as 
compared to the critical value 5c [l8l |: 
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• (5i < Sc- In this case an event horizon always exists 
on each side of the wormhole for all values of Ng. 
We could say that these solutions behave somewhat 
like Schwarzschild black holes. 

• (5i > Sc'. The structure of horizons is more com¬ 
plicated and (on each side of the wormhole) one 
can find two, one (degenerate), or no horizons, like 
in the usual RN solution of GR. Nevertheless let 
us stress that in all these cases the structure close 
to the center undergoes important changes as com¬ 
pared to their GR counterparts. 

• (5i = Sc'. If JVg > Nc, one finds that there are 
two horizons located symmetrically on each side 
of the wormhole. If Ng = Nc, the two horizons 
meet at the wormhole throat, r = Vc (or a; = 0). 
If Ng < Nc then the horizons disappear yielding 
a kind of black hole remnant. The existence of 
such remnants, which can be originated as the end 
state of a black hole under Hawking evaporation or 
due to large density fluctuations in the early Uni¬ 
verse [13 , might be of special relevance for the un¬ 
derstanding of the information loss problem [35l| . 
Besides, they have potential observational conse¬ 
quences [13 • 


III. EUCLIDEAN EMBEDDINGS 

The divergence of the metric function (1141) as r —>■ Tc 
when ^ 5c also implies curvature divergences there. 
However, the wormhole structure and physical properties 
such as total charge, mass, and density of field lines are 
finite and as well-behaved as in the case di = 5c, which is 
completely free from curvature divergence^ This sug¬ 
gests that curvature divergences might not be as trou¬ 
blesome as they seem to be in structureless scenarios, 
i.e., when they occur at a point rather than around a 
finite-size topological structure such as a wormhole. To 
get an intuitive idea of the differences and similarities 
between the smooth case = 5c and the divergent case 
(5i ^ 5c, we find it useful to construct an Euclidean em¬ 
bedding of the spatial equatorial sections of these geome¬ 
tries. This can be done by considering the 0 = 7r/2 and 
t =constant section of the line element (I5|) expressed in 
terms of dx^ = dr^ /cr_, which yields 

dl^ = -h V , (15) 


^ We note that in GR electrovacuum scenarios resulting from non¬ 
linear theories of electrodynamics, such as Born-Infeld theory 
[^ . the metric may be finite at r = 0 but nevertheless have 
curvature divergences at that point [s^l . As already mentioned, 
singularity avoidance in such a context is done at the price of 
violations of the energy conditions and/or ill-definiteness of the 
underlying electromagnetic theory [^. 


and embedding it into a three-dimensional Euclidean 
space with cylindrical symmetry of the form 


dl^ + (16) 


One just needs to find the function ^(r) that leads to the 
line element (fT5]l . Since the region relevant to our discus¬ 
sion is the neighborhood of r ^ rc, where the wormhole 
throat is located, we can take the near wormhole expan¬ 
sion da together with a- ~ 4(r — Tc)lrc to get 



To illustrate this procedure, we will restrict ourselves to 
the simplest cases, namely, i) regular, horizonless black 
hole remnants (corresponding to (5i = 5c and Ng < Nc) 
and ii) the Reissner-Nordstrdm-like case, > 5c, in those 
cases in which A{x) > 0 near the wormhole. This lat¬ 
ter case represents both configurations without horizons 
{naked) and configurations with two horizons. In general, 
for a line element of the form dP = C{r)dr^ +r^dVP, the 
function ^(r) must be of the form = C{r) — 1. Since 
in our case the functions C{r) both diverge as r —>■ r^, 
we can approximate ~ C'(r). We thus get that 




I (^c-A^G 
^ iNa 


^c\/r - Tc 


if = 5c 


(5i-5c) = 



if > 5c 


(18) 


In Figs. [2] and [3] we have plotted the resulting embed¬ 
dings for horizonless solutions. In the case of Reissner- 
Nordstrom-like solutions with two horizons a representa¬ 
tion is also possible, but the inner horizon and the worm- 
hole surface are very close to each other, which would re¬ 
quire a distortion of the radial coordinate for its graphical 
representation. In both cases the presence of the worm- 
hole structure becomes manifest. We note that, while in 
the regular = 5c case the region around the worm- 
hole throat is completely smooth, for the > 5c case 
it presents a cusp. See also Fig. |3]for a 2—dimensional 
comparison. 

To see in a more quantitative way the differences be¬ 
tween those conhgurations, we consider the Kretschmann 
scalar corresponding to the Euclidean surfaces repre¬ 
sented in Figs. [2] and [3l which can be computed with 
the line elements given in (1181) using the formula K 2 D = 

72 ^ ■ The result is 


K 2 D = 


64(N^-N„f 1 


if (5i = 5c 


(Si-Sa f 1 

N'^ 4S'fS^ rc(r—rc)r‘^ 


if (5i > 5c 


(19) 


This puts forward that two apparently similar surfaces 
can have very different properties as far as curvature 
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Figure 2. Euclidean embedding of the 6 = n/2 spatial section 
of a regular wormhole (5i = 5c)- The vertical axis represents 
the function ^(r). 


Figure 3. Euclidean embedding of the 0 = 7r/2 spatial section 
of a wormhole with curvature divergences at its throat {naked 
singularity with 5i > Sc)- The vertical axis represents the 
function ^(r). 


scalars are concerned. From the definition of K 2 D one 
readily finds that for C{r) = {r — rc)°', the geometry is 
free from curvature divergences if a < — 1 and diverges 
otherwise. The regular configuration (5i = Sc saturates 
the bound a = — 1. Similarly, one can also verify [Tsj l 
the finiteness of other curvature invariants for the regu¬ 
lar cases. 


IV. CONFORMAL DIAGRAMS 

We have already discussed the horizons of the geome¬ 
try from Eg. (1141) . In order to draw a conformal diagram 
of the geometry we have to look into the nature of the hy¬ 
persurface a: = 0 (or r = re), where the wormhole throat 
is located. A hypersurface S' in a manifold M is called 
space-like, null, or time-like when the tangent space to S 
at each point has this same character. Therefore, the 
normal to a space-like, time-like, or null hypersurface 
must be time-like, space-like, or null, respectively. For 
a given surface defined by f{v,x,9,(j)) = /o = constant, 
the normal vector is defined as = ^r- For the hy¬ 
persurfaces X = /o, the normal is = (0,1,0,0), which 
yields = Acr^. Therefore, if A > 0 the hypersur¬ 

face is time-like, if A < 0 it is space-like, and if A = 0 
then it is null. Accordingly, looking at Ea. dHl) we can 
distinguish the following cases: 

• (i5i < (5c) => a; = 0 is a space-like hypersurface. 





Figure 4. Representation of ^(r) as a function of r. The 
dashed curve represents the regular wormhole configuration 
(5i = Sc), while the continuous curve is one of the solutions 
with curvature divergences. The curves have been normalized 
to make them coincide at 2 = 10. Note that both curves are 
continuous everywhere and have divergent derivative at 2 = 1. 


• (5i > (5c) => a: = 0 is a time-like hypersurface. 

• (5i = (5c, Nq > Nc) => x = 0 is a space-like hyper¬ 
surface. 

• (5i = dc, Nq = Nc) => X = 0 is a null hypersurface. 

• (5i = 6c, Nq < Nc) X = 0 is a time-like hyper¬ 
surface. 
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Note that this classification of the surface x = Q (or 
r = Tc) differs from that initially given in [l 8 j |. where 
the metric was represented using the function r as a co¬ 
ordinate. As pointed out above following Eq.®, r is 
not a valid coordinate at r = rc, which invalidates the 
classification and some aspects of the Penrose diagrams 
provided in [T^. 

Taking into account the number and type of horizons, 
together with the existence or not of curvature diver¬ 
gences, one finds seven different possible causal struc¬ 
tures. In this sense, regular configurations without hori¬ 
zons appear in Fig. [5l with horizons in Fig. El and with 
a horizon coinciding with the wormhole throat in Fig. 0 
In Fig. [ 8 ] we have represented the Schwarzschild-like so¬ 
lutions (a single non-degenerate horizon) and Figs. I^ITUl 
and EH are the Reissner-Nordstrom-like solutions corre¬ 
sponding to naked singularities, extreme black holes, and 
two-horizon black holes, respectively. 



Figure 7. Penrose diagram for the case 5i = Sc, Nq = Nc- In 
this case the event horizons coincide with the location of the 
wormhole throat, which becomes a null hypersurface. 



Figure 5. Penrose diagram for the case = Sc, Nq < Nc- 
The wormhole follows a time-like trajectory. 



Figure 6. Penrose diagram for the case = Sc, Nq > Nc. 
The wormhole represents a space-like hypersurface. 



Figure 8. Penrose diagram for the case < Sc- The wormhole 
is a space-like hypersurface with curvature divergences at r = 

Tc- 

large values of the affine parameter (geodesic complete¬ 
ness). Comparison of some relevant results with their GR 
counterpart (the standard Reissner-Nordstrom solution) 
is provided. 

A. Geodesics in a metric-afHne spacetime 

Geodesics are curves whose tangent vector is paral¬ 
lel transported along itself. This definition generalizes 
the concept of “straight lines” of Euclidean geometry to 
curved geometry. A geodesic curve 7 ^ = x^{\) with tan¬ 
gent vector and affine parameter A satisfies, in 

a coordinate basis, the following equation @ 


V. GEODESICS 

In this section we present a detailed description of 
the geodesic behaviour corresponding to different con¬ 
figurations of our wormhole solutions, depending on the 
charge-to-mass ratio and the number of charges involved. 
Our goal is to determine whether geodesic curves cross¬ 
ing through the wormhole can be extended to arbitrarily 


or, equivalently 

dx“ Ax^ 
dA2 dA dA 


( 20 ) 


( 21 ) 


which is a set of second-order differential equations. 
These equations have a unique solution for a given con¬ 
nection T))^ and initial conditions a:^( 0 ), (dx^/dA)|o. 









Figure 9. Penrose diagram for the case > Sc (in GR this 
represents a naked singularity). 



Fignre 10. Penrose diagram for the case tfi > Sc (one extremal 
horizon). 


In metric-affine theories like the ones considered here, 
one assumes the a priori existence of independent metric 
and affine structures. As is wefl-known, the metric (or 
causal) structure can be used to define an affine structure 
in terms of the Christoffel symbols of the metric (Levi- 
Civita connection). This determines a set of geodesics 
which, according to the Einstein Equivalence Principle 
(EEP), coincide with the paths followed by test parti¬ 
cles. These curves also extremize the length between 
its endpoints In the metric-affine case, the indepen¬ 
dent connection can also be used to define a different 
set of geodesic paths. If the theory is constructed as¬ 
suming the EEP, i.e., not coupling the connection to the 
matter fields, then the independent connection should 
just be a gravitational field which contributes to gen¬ 
erate the spacetime metric but is expected not to act 
directly on the matter fields [d^. If, on the contrary, 
the matter fields are coupled directly to the indepen¬ 
dent connection in the action, one should then study 



Figure 11. Penrose diagram for the case (5i > Sc (two hori¬ 
zons). 

those geodesics as physically meaningful. In our case, 
the geometry has been derived assuming the existence 
of just an electric field, which is insensitive to the de¬ 
tails of the (symmetric) connection. In fact, variation 
of the Maxwell action leads to = 0, and 

given that V = d^^—g — this equation boils 
down to i'/^F^’') = 0, which has no dependence on 
the particular connection P^,^ used to define the covari¬ 
ant derivative. For this reason, our focus will be on the 
geodesics of 

Instead of considering the geodesic equation itself to 
obtain the paths followed by test particles, it is more 
convenient to exploit the symmetries of the problem to 
obtain conserved quantities that simplify the analysis 
S13 . To proceed, we note that the geodesic equations 
(l?n) can be derived from an action principle S = J dXL 
with Lagrangian L = \gnv^^^§\ ■ The momenta asso¬ 
ciated to this Lagrangian are given by = dL/dx^, 
with which leads to ■ One thus 

finds that H = Pf^x^ — L = L, which in terms of the 
momenta reads as H = ip^'^(x)p^pi,. Note that the ab¬ 
sence of a potential term in this Hamiltonian puts for¬ 
ward that geodesic trajectories can be seen as represent¬ 
ing free particles in a curved geometry. The Hamiltonian 
equations of motion are thus x^ = dH/dp^ = g^'^Pu and 
p^ = —dH/dx^ = —^PaPgd^g^^^. Using these equa¬ 
tions, one can easily verify that dx^/dX reproduces the 
geodesic equation (1^ . It is also a trivial matter to show 


We note that in Palatini theories the situation is not as simple 
as generally thought because even if one assumes the postulates 
of metric theories of gravity in the construction of the theory, 
violations of the EEP are still possible 0- 
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that dH/dX = 0, which implies that the Hamiltonian 
is a conserved quantity. Given that does not de¬ 
pend explicitly on the coordinates t and ip, one finds 
that pt = 0 and = 0 represent other two conserved 
quantities. In terms of the line element (HI), we thus have 
that dt/dX = E/A and dip/dX = L/r^, with E and L 
constants, where we have taken 9 = tt/2 because due to 
spherical symmetry the geodesics must lie on a plane. 
If one uses (HI then E = For timelike 

geodesics, E can be interpreted as the total energy per 
unit mass, and L as the angular momentum per unit 
mass. For null geodesics E and L lack meaning by them¬ 
selves, since it is not possible to normalize the tangent 
vector, but the quotient L/E can be interpreted as the 
apparent impact parameter in the asymptotically fiat in¬ 
finity. By rescaling the affine parameter by a constant, 
the Hamiltonian can be set to =Fl for time-like/space-like 
geodesics, respectively, whereas for null geodesics E[ = 0. 
The constancy of the Hamiltonian, therefore, allows us 
to write the following constraint for the geodesic tangent 
vector: 



where /c = 0 for null geodesics and k = I for time¬ 
like geodesics. For time-like geodesics, A represents 
the proper time of the particle following the geodesic, 
whereas for null geodesics it is an affine parameter. Us¬ 
ing the conservation relations, (|2^ turns into 


1 


a 


2 

+ 




(23) 


Under a rescaling of the form dy = dx/a+, (1^ can be 
seen as a single differential equation akin to that of a 
classical particle in a one dimensional potential of the 
form 


V{x) = A 



(24) 


Had we used the line element o, which is also valid 
in case of having event horizons, Eq.([23| would still be 
valid. From now on, we will study the behaviour of the 
geodesics in terms of the potential. Since V(x) is a func¬ 
tion of r(x), which is even in the variable x, it turns out 
to be also an even function. Our description of the po¬ 
tential will thus be restricted to the a; > 0 sector, which 
has a direct correspondence with the GR case. 


B. Radial null geodesics 


Radial null geodesics are characterized by k = 0 and 
L = 0 and, therefore, satisfy the equation 


1 


(T 


2 

+ 



= E^ 


(25) 


which is insensitive to the details of the function A. Using 
the relation between r and x, one can write x^ = r^CT_, 
which implies dx/dr = ±ct+/ct_ , with the minus sign 
corresponding to x < 0. This turns (P5| into 


I 

cr_ 



(26) 


This last equation admits an exact solution of the form 

[ 2Fi[-i,i,f;^]r if X > 0 

±E-X{x) = l 

[ 2x0 - 2Fi[-i, i, ^]r if a; < 0 

(27) 

where 2 F 1 [a, b, c; y] is a hypergeometric function, xq = 

~ 0.59907, and the ± sign 
corresponds to outgoing/ingoing null rays in the x > 0 
region. It should be noted that given that dr/dX is a 
continuous function, the solution (l?7| is also unique. For 
X —> 00 the series expansion of this expression yields 
EX{x) Ki r + 0{r~^) ~ x and naturally recovers the GR 
behavior for large radii (see Fig fT^ . As the wormhole 
throat is approached, one finds EX(x) ~ xq ± \/r — Xc ~ 
xo-l-x/2, with the -|- (—) sign corresponding to the branch 
with X > 0 (x < 0). Numerically one verifies that this 
approximation is very good within the interval x/r^ G 
] — 1,1[. In the limit x —>■ —00, A(x) « x -|- 2xo recovers 
the linear behavior but shifted by a constant factor. 


A 



Figure 12. Affine parameter X{x) as a function of the radial 
coordinate x for radial null geodesics (outgoing in x > 0). In 
the GR case (green dashed curve in the upper right quadrant), 
A = X is only defined for x > 0. For radial null geodesics in 
our wormhole spacetime (solid red curve), A(x) interpolates 
between the GR prediction and a shifted straight line A(x) « 
X -I- 2xo, with Xq « 0.59907. In this plot E = 1 and the 
horizontal axis is measured in units of Xc. 

It is remarkable that the affine parameter A = A(x) 
given in dSZI extends over the whole real axis. This 
contrasts with the GR prediction for electrovacuum con¬ 
figurations, where null radial geodesics take the form 
(dr/dA)^ = E^ and whose solution for outgoing/ingoing 
geodesics is of the form r(A) = ±EX. In the GR 
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case, therefore, the affine parameter A is only defined on 
the positive/negative (outgoing/ingoing) side of the real 
axis because the function r(A) is positive definite. The 
Schwarzschild and Reissner-Nordstrom black holes in GR 
are thus said to be geodesically incomplete as far as null 
geodesics are concerned. In our case, on the contrary, 
radial null geodesics are complete and this occurs for ar¬ 
bitrary choices of the parameter Ji. This is relevant be¬ 
cause generically a curvature divergence occurs at a; = 0, 
where the wormhole throat is located. Onl y fo r the case 
i5i = Sc is the geometry completely regular [l^. Ea. (l27l) . 
therefore, puts forward that radial null geodesics are the 
same for all the wormhole configurations, regardless of 
the possible existence of curvature divergences. We also 
note that the radial null geodesics of the metric are 
the same as those corresponding to the auxiliary metric 
defined in dm- 


vergent barrier as a; —>■ 0 (see FiglTS] case A for 
details). Note that, unlike in the case of a particle 
in a potential, in our case the parameter > 0 
always. This means that no stable photon orbits 
may exist at the minimum of the potential in the 
dotted curve of plot A in FiglTSl which occurs in a 
region where V(x) < 0. 

• If = Sc (regular case), the potential is regular 
at a; = 0 (see FiglT^ case B for details), behaving 
there as 


V{x)/L^ 



(29) 


C. Null geodesics with L ^ 0 

For null geodesics (k = 0) with angular momentum 
L 7 ^ 0 , a geodesic coming from r —>■ oo (or, equivalently, 
X — >■ ±00 ) starts seeing the typical centrifugal barrier 
term of GR, which grows from zero like V m This 
barrier is positive and negligible far away but grows as 
the center is approached. The behavior as a; —>■ 0 depends 
crucially on the parameters (5i and S 2 = SiNc/Ng that 
characterize the background geometry. In fact, in the 
limit a; —>■ 0 , we have 

V{x)^-^^-b ( 28 ) 

with a = («:± fand b = (« ± §) 

This leads to the following cases when k = 0: 

• If > Sc (Reissner-Nordstrom-like case) one finds 
an infinite repulsive barrier as x = 0 which makes 
all geodesics bounce at some r > Xc, preventing 
them from reaching the wormhole in much the same 
way as it happens in the usual Reissner-Nordstrom 
solution of GR Q , where L 7 ^ 0 null geodesics can¬ 
not reach the central singularity. For certain values 
of the charge, the potential may have a local max¬ 
imum and then a minimum before reaching the di- 


(recall that S 2 = SiNc/Nq). The potential has an 
extremum at x = 0, which can be a minimum in 
between two local maxima for some values of the 
parameters. All geodesics with energy greater than 
the maximum of the potential will go through the 
wormhole (see FiglM] for details). Stable photon 
orbits are possible at the minimum of the potential 
if Nq < Nc because then V (0) > 0. Bounded orbits 
can also exist near the wormhole if a photon is emit¬ 
ted with 0 < < Vmax, being Vmax the maximum 

value of the potential (due mainly to the centrifu¬ 
gal barrier). A glance at Fig ITS] shows that bounded 
photon orbits can exist around the wormhole throat 
if a photon is emitted inside this region (not coming 
from infinity) with an E^ smaller than the potential 
barrier. The dotted potential in FiglTSl indicates 
that such photon trajectories would be bounded 
by the centrifugal barrier. These geodesics can get 
into a black hole region after crossing an external 
event horizon, go through the wormhole, and get 
out of the black hole region after crossing the other 
event horizon (recall that, from the definition of 
V(x), the zeros of V(x) coincide with the zeros of 
A(x), which signals the presence of horizons). This 
photon would then bounce due to the centrifugal 
barrier and enter the black hole region to repeat 
the process in reversed directioi{l. 


• If (5i < Sc (Schwarzschild-like case) the potential 
becomes infinitely attractive at x = 0, with the pos¬ 
sibility of having a maximum before that point, de¬ 
pending on the number of charges Nq. All geodesics 


with energy above that maximum hit the wormhole 
(see FiglT3]case C for details). With the approx¬ 
imate form of the potential in the x —> 0 region, 
namely. Eg. (1^51) . one can verify that 


dA 



(b + E^) 
4a i 


® It might be useful at this point to recall the discussion about 
the advanced Eddington-Finkelstein coordinate provided in the 
paragraph containing Eq. ®. 


dx 20^/2 


(30) 
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Figure 13. Representation of the effective potential for null geodesics with L = 1. Plots A, B, and C correspond to the charge- 
to-mass ratios = 1.5(5c, di = dc and = 0.3Sc, respectively, where the curves represent the cases with Ng = l,Ac,8A^c 
(solid curve, dashed curve, and dotted curve, respectively). Plots D, E, and F provide a comparison of different values of 
for Ng = 1 (plot D), Ng = JVc (plot E), and JVg = 8Nc (plot F), being = 1.55c the dashed (red) curve, 5i = 0.35c the dotted 
(blue) curve, and 5i = 5c the solid (green) curve. The same colors have been used in plots A, B, and C to represent the value 
of 5i. 


which leads to 


A(a;) 


X X i 
3 a 


(. Hb + E^) 

V 10 



(31) 


where the integration constant has been chosen so 
as to make A(0) = 0. Note that despite the diver¬ 
gence of the potential at a; —>■ 0, the affine parame¬ 
ter is smooth across that point, as shown in FiglThl 
In fact, given that the right-hand side of dX/dx 
is a smooth function, the solution (ED) turns out 
to be unique once initial conditions are specified. 
This confirms that null geodesics with L ^ 0 are 
also complete in this spacetime. Note in this sense 
that L ^ 0 null geodesics in Schwarzschild space- 
time are not complete because r = 0 is reached 
in a finite affine time and there is no way to ex¬ 
tend the affine parameter to an hypothetical region 
r < 0. Another way to see the incompleteness of 
these geodesics in GR is through the conservation 
of angular momentum equation, L = which 

makes the angle tp undefined as r —>■ 0. In the 
wormhole case, on the contrary, the angular veloc¬ 
ity is finite at the wormhole throat r = Tc, which 
avoids this problem too. Similarly as in the 5i = Sc 
case with Nq > Nc, bounded photon trajectories 
with 0 < < Vmax can exist in the region close 

to the wormhole. 





Figure 14. Maximum of the potential for null geodesics with 
5i = 5t:, against the number of charges. All photons emitted 
from infinity with E/L > 0.5 will be able to go through the 
wormhole. 


D. Radial time-like geodesics 

In the radial time-like case {k = 1, L = 0), the be¬ 
havior far away from the wormhole (on both sides) is 
identical to that found in GR, being dominated by an 
attractive potential V{x) ~ 1 — As x —>■ 0, the 
potential is dominated by the approximation (1281) . The 
dependence on 5i, therefore, determines the evolution of 
the geodesics: 

• In the RN-like case, 5i > 5c, there is an infinite 
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that there exist bound orbits for particles with 0 < 
< 1, having a stable point near the wormhole if 
^min ^ b- 

• For i5i = Sc the potential is finite at a; = 0, having 
the form 




+ 4'“- 


2Ng 

3Nc 


7 N, 
240 7V< 


Q 4 
X 


(32) 


Figure 15. Representation of the potential V{x) for null 
geodesics in the regular case = Sc- The solid curve rep¬ 
resents the case Nq = 1, the dashed curve corresponds to 
Nq = Nc, and the dotted curve (with two event horizons lo¬ 
cated at the zeros of V{x)) represents the case Nq = 1.2Nc- 



Figure 16. Representation of the affine parameter \(x) for 
geodesics with = 0.35c, and Nq = 2Nc- The solid (red) 
curve represents the approximation iu, the dashed (green) 
curve represents the null case with L = 1, and the dotted 
(blue) curve is the time-like case with L = 0. The potential 
that generates the dashed curve is similar to the dashed curve 
of plot C in Fig |13l The potential of the dotted curve is similar 
to the dotted potential of plot C in Fig |17l 

repulsive barrier at a; = 0, which arises right after a 
minimum. Massive particles with L = 0, therefore, 
cannot reach the wormhole (this is similar to GR, 
where massive particles cannot reach the central 
singularity 0)- From plot A of FiglTTl one finds 


If Nq < ^Nc, then a: = 0 is a minimum of the po¬ 
tential. This means that massive particles can stay 
at rest at a; = 0 if < Nc, because then R( 0 ) > 0 . 
For larger values of the charge, the wormhole throat 
a; = 0 is hidden behind an event horizon (one on 
each side) and 14(0) < 0, which prevents the exis¬ 
tence of stationary points. Bound orbits may exist 
in this region as long as E'^ is smaller than the 
maximum of the centrifugal barrier, similarly as in 
the case of photons with L yf 0. Therefore, a mas¬ 
sive particle whose E'^ drops below the maximum of 
the potential in the region near the wormhole can 
oscillate around the wormhole bounded by the cen¬ 
trifugal barrier. This oscillatory motion would be 
possible even when Nq > Nc, i.e., when the worm- 
hole is hidden by event horizons (one on each side). 
Note that for Nq > Nc we cannot have the parti¬ 
cle at rest at the minimum of the potential because 
14(0) <0 can never satisfy dx/dX = 0. This is con¬ 
sistent with the fact that a massive particle within 
the event horizon cannot stay at rest. 

• In the Schwarzschild-like case, i5i < Sc, there is 
an infinite attractive well as a; —>■ 0. All radial 
time-like geodesics, therefore, reach the wormhole. 
From the approximate form of the potential in this 
region, one can verify that (I5T]) provides a good de¬ 
scription of the geodesics around x = 0. Therefore, 
the affine parameter can be smoothly extended 
across x = 0 also in the time-like case (see Fig. 
M for a comparison of this case with the null and 
approximated cases). Bounded orbits around the 
wormhole also exist in this case. 


E. Time-like geodesics with L 0 

For nonzero angular momentum time-like geodesics 
obey a potential which is the sum of the two previous 
cases. The qualitative features of the previous cases are 
also manifest here when L > 0. For (5i < 5c the approx¬ 


imate formulas developed in the null case for the affine 
parameter near x = 0 are also valid. In fig. [TS]we have 
plotted the potentials in the particular case L = 2 in the 
three cases 5i > 5c, 5i = Sc and 5i < 5c and different val¬ 
ues of the number of charges Nc- We can thus conclude 
that time-like geodesics are also complete in our worm- 
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D: L=0 and Nq=1 


E: L=0 and Nq^N^ 


F: L=0 and N,=2Nq 


Figure 17. Representation of the effective potential for time-like radial geodesics (with L = 0). In plots A, B, and C the curves 
represent cases with different values of the charge parameter Nq (the solid curve is Nq = 1, the dashed curve has > 1, and 
the dotted curve has the largest charge). Plots D, E, and F provide a comparison of different values of for A^q = 1 (plot D), 
Nq = Nc (plot E), and Nq = 2Nc (plot F), being Si = 1.55c the dashed (red) curve, 5i = 0.35c the dotted (blue) curve, and 
5i = 5c the solid (green) curve. The same colors have been used in plots A,B, and C to represent the value of 5i. 


hole spacetime, which clearly contrasts with the results 

I 


F. Stationary null orbits 

A stationary orbit occurs when the energy of a particle 
coincides with the potential energy at an extremum of 
the potential, i.e., when = V{xij) and Vx\x=xo = 0. If 
this extremum is a minimum, a slight perturbation will 
make the orbit oscillate around the minimum. If it is a 
maximum, then the stationary point is unstable. 

Since for given values of the parameters 5i and 62 our 
geometry for \x\ 1 is locally indistinguishable from 

that provided by GR, the stationary orbits that one finds 
at large radii are almost coincident with those found in 

I 


VI. SUMMARY AND CONCLUSIONS 

In this work we have studied geometrical aspects of 
a family of wormhole solutions supported by a spheri¬ 
cally symmetric electric field which arise in high-energy 
extensions of Einstein’s theory formulated a la Palatini. 


for Schwarzschild and Reissner-Nordstrom black holes of 
GR. 


the standard Reissner-Nordstrom solution (see FigHH). 
The stationary orbits that occur near the wormhole de¬ 
part from those found in GR as we get closer to a; = 0. In 
Figs. EDI EH and E2] we plot the location of the station¬ 
ary null orbits as a function of the number of charges for 
both the Reissner-Nordstrom solution of GR and for our 
wormhole geometry and for various values of 5i, equal 
to, lower and greater than 5c, respectively. One can ver¬ 
ify that for certain values of the charge, the stationary 
orbits may not exist in the GR case but persist in the 
wormhole scenario. Note that stationary null orbits exist 
at the wormhole throat for Nq < Nc when 5i = 5c (Fig. 

EQl). 


Euclidean embeddings have been used to illustrate that 
similar wormhole structures may possess very different 
properties as far as curvature scalars are concerned. Gon- 
formal diagrams of the different wormhole configurations 
have been provided to update preliminary analyses car¬ 
ried out in . It should be noted that the conformal di- 
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Figure 18. Representation of the effective potential for time-like geodesics with L — 2. In plots A,B, and C the curves represent 
cases with different values of the charge parameter Nq (the solid curve is Nq = 1, the dashed curve has > 1, and the dotted 
curve has the largest charge). Plots D, E, and F provide a comparison of different values of 5i for Nq — 1 (plot D), Nq = Nc 
(plot E), and Nq = 4:Nc (plot F), being = 1.5(5c the dashed (red) curve, = 0.3(5c the dotted (blue) curve, and Si — Sc the 
solid (green) curve. The same colors have been used in plots A, B, and C to represent the value of di. 



B: di=dc. L=1/2,1,2 and Nq= 2Ne 


V(x) 



Figure 19. Comparison of the effective potential for null geodesics in the GR and wormhole cases. The solid (blue) curves 
represent the wormhole case and the dashed (orange) curves the GR case. The GR potential is only defined for x > 0. Note that 
the convergence between the GR and wormhole cases occurs very quickly for not very large values of x > 0, which manifests 
the non-perturbative nature of the wormhole geometry. Plot A represents the Reissner-Nordstrom like case (di > 5c). Plot B 
corresponds to the regular configurations (5i = 5c). Plot C is the Schwarzschild-like case (5i < 5c). 


agrams containing curvature divergences can be extended 
to include the region across the wormhole. By doing this, 
Figiniwould look like Figl5]with the straight diagonal line 
replaced by a zig zag line. Similar modifications would 
be necessary in FigsfTUlandITT] 

We have carried out a detailed study of the geodesic 
structure of these spacetimes finding that the three 
possible configurations, namely, Reissner-Nordstrom-like 
(5i > 5c), Schwarzschild-like (5i < 5c), and Minkowski- 
like (5i = 5c) are geodesically complete. This is so despite 
the fact that only in the case 5i = 5c are curvature scalars 
regular everywhere. In the other cases, 5i ^ Sc, curvature 


divergences appear at the wormhole throat. This result 
puts forward, through an explicit example, that the blow 
up of curvature invariants such as the squared Ricci ten¬ 
sor or the Kretschmann scalar does not necessarily imply 
geodesic incompleteness, which is the principal criterion 
to determine if a spacetime is singular or not [J ■ We thus 
conclude that the family of geonic wormhole solutions 
(with or without event horizons) provided by the Pala¬ 
tini version of quadratic gravity and/or the Born-Infeld 
theory of gravity represent non-singular spacetimes. Re¬ 
markably, this result follows from the interplay between 
the Palatini gravity model and the standard Maxwell 
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Figure 20. Radius of the stationary null orbits in units of 
Tc (vertical axis) against the number of charges (horizontal 
axis) for a black hole with di = Sc. The blue (upper) and 
green (flat) solid lines are for the wormhole configuration, the 
dashed (orange) line is for the RN black hole of GR. Notice 
that the stable (flat) branch of stationary orbits ends at Nq = 
Nc- The upper solid (blue) curve smoothly tends to the GR 
prediction for large values of Nq. 



Figure 21. Radius of the stationary null orbits in units of Vc 
(vertical axis) against the number of charges (horizontal axis) 
for a black hole with = 0.5 Sc. The solid (blue) line is for 
the wormhole, the dashed (orange) line is for the RN black 
hole of GR. 


field, not from the introduction of exotic energy sources 
in the framework of GR [l^ . 

We would like to stress that the wormhole (topological) 
structure of our spherically symmetric solutions is the 
crucial element that avoids geodesic incompleteness [d^ . 
The case of radial null geodesics in the Schwarzschild- 
like case (di < Sc) is very illustrative to understand this 
point. For ingoing null or time-like geodesics, as the 
time V in m passes by, the radial coordinate x must 
decrease [see the discussion around ®]- This coordinate 
goes from -|-oo to —oo while the radial function r(x) al¬ 
ways remains positive. A spherical shell of particles or 
radiation going in from x —>■ -boo is seen to collapse into 
a minimal surface of area A = 47rr^ at a; = 0 before 
bouncing off as an outgoing shell of particles/radiation 
as the wormhole is crossed in the direction of x —>■ — oo. 
In the case of GR, the same shell of particles/radiation 
would have reached the center r = 0 in a finite affine 
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Figure 22. Radius of the stationary null orbits in units of 
(vertical axis) against the number of charges (horizontal axis) 
for a black hole with (5i = 1.05 * ic The solid blue and green 
lines are for the Palatini black hole, the dashed (orange) line 
is for the RN black hole of GR. 


parameter with no possible extension beyond that point 
(because r is always positive, r = 0 represents its min¬ 
imum value, and there is no possibility to go back to 
larger values of r within the event horizon). If angular 
momentum is considered, the situation worsens in GR, 
as the angular velocity dip/dX = Ljr'^ diverges as r —0 
and the hypothetical extensions beyond that point would 
have completely undetermined the angular coordinate ip. 
In the wormhole case, (p is well defined at r = Tc, which 
guarantees its smooth continuation across a; = 0. The 
Euclidean embeddings of Sec lIIIl can be used to visualize 
how any smooth curve (such as spatial geodesics satis¬ 
fying dip/dX = L/r^) reaching the wormhole throat can 
be continued to the other side despite the possibility of 
having curvature divergences at the throat. 

Before concluding, we note that our analysis has fo¬ 
cused on the properties of individual geodesics. Since 
physical observers are sometimes represented as congru¬ 
ences of geodesics, there remains to determine how con¬ 
gruences behave as they approach regions with curvature 
divergences. This point, in fact, has been used in the lit¬ 
erature to classify the strength of curvature singularities, 
by considering the behaviour of the volume element asso¬ 
ciated to a physical observer travelling through the sin¬ 
gularity, to determine whether it is crushed/ripped apart 
in the process or can safely cross it [d^ - lsoj] . In addition, 
one could also consider the fundamental wave-like nature 
of particles and test the singularity by quantum scatter¬ 
ing experiments [5l| . A detailed study of these aspects is 
currently underway and a preliminary analysis has been 
reported in [5^. 
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